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Abstract 

With heavy quark limit and hierarchy approximation Xqcd ^ m r> "C ito.b, we analyze the 
B — > D°D° and £? s — > D°D° decays, which occur purely via annihilation type diagrams. As 

a roughly estimation, we calculate their branching ratios and CP asymmetries in Pcrturbative 

D ' n — ° _« 

QCD approach. The branching ratio of B — > Z) U D is about 3.8 x 10 that is just below the 

latest experimental upper limit. The branching ratio of B s — * Z? Z?° is about 6.8 x 10~ 4 , which 



o 



X 



could be measured in LHC-b. From the calculation, it could be found that this branching ratio is 

(N 
oo 

because of interference between weak and strong interaction. However, these asymmetries are 
too small to be measured easily. 



1 Introduction 

In the Standard Model (SM), CP-violation (CPV) arises from a complex phase in the Cabibbo- 
Kobayashi- Maskawa (CKM) quark mixing matrix, and the angles of unitary triangle are defined 
as [I]: 



(3 = arg 



VgVu 



a = an 



V t tV td 



7 = arg 



V ub V ud 

viv cd 



1) 



v* b v ud 

In order to test SM and search for new physics, many measurements of CP-violation observables 
can be used to constrain these above angles. It is well known that we measure (3 precisely using 
the golden decay mode B — > J/ipK s ; the angle a can be determined with decay B — > hit and 7 
could be measured precisely in Large Hadron Collider (LHC) with decay mode B s — > D S K. 



* e-mail: liying@ytu.edu.cn 



1 




(a) (b) 



Figure 1: The quark level Feynman diagrams for B& —* D°D° process 



Besides the above channels mentioned, many other channels are used to cross check the mea- 
surements. Among these decays, B — > DD decay is considered to test the j3 measurement. For 
B — > DD decay, the analysis based on SU(3) symmetry [2], iso-spin symmetry |3j and factorization 
approach [1] have been done in last several years. However, the the calculation of decay B° — > D°D° 
has difficulties. It is a pure- annihilation diagram decay, also named W-exchange diagram decay, 
which is power suppressed in factorization language. The quark diagrams of this decay are shown 
in Figure [TJ Theoretically, QCD factorization approach (QCDF) [5] and soft collinear effective 
theory (SCET)[6] can not deal decays with two heavy charmed mesons effectively. In Ref.jjTl [5J, 
perturbative QCD (PQCD) has been exploited to B meson decays with one charm meson in the 
final states and the results agree with experimental data well. Specially, the pure annihilation type 
B decays with charmed meson were studied in Ref.[8j. 

In the standard model picture, the W boson exchange causes bd — > cc, and the uu quarks are 
produced from a gluon. This gluon attaches to any one of the quarks participating in the W boson 
exchange. In decay B — > D°D° , the momentum of the final state D meson is ^mg(l — 2r 2 ), with 
r = mo/niB- If we consider heavy quark limit and hierarchy approximation Xqcd *C mo <C tob, 
the D meson momentum is nearly mg/2. According to the distribution amplitude used in Ref.[7], 
the light quark in D meson carrying nearly 40% of the D meson momentum. So, this light quark 
is still a collinear quark with 1 GeV energy, like that in B -> DM IS], B -> K(ir)ir [TO] 
decays. The gluon could be viewed as a hard gluon approximatively, so we can treat the process 
perturbatively where the four-quark operator exchanges a hard gluon with uu quark pair. Of 
course, we are able to calculate the diagrams if charm quark and up quark exchange. As a roughly 
estimation, we give the branching ratio and CP-violation of B^ s ~^ D°D° . 

In this article, the analytic formulas for the decay amplitudes will be shown in the next section. 
In section [3j we give the numerical results and summarize this article in section HI 
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2 Analytic formulas 

For simplicity, we set B meson at rest in our calculation. In light-cone coordinates, the momentum 
of B, D° audi? are: 

Pb = ^§(1, 1, 0); P 2 = ^§(1 - r\ r\ 0); P 3 = ^§(r 2 , 1 - r\ 0). (2) 

we define the light (anti-)quark momenta in B, D° and D° mesons as k\, k 2 , and k 3 as: 

h = ( Xl P+,0,k 1T ), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 P 3 ~,k3r). (3) 

In PQCD, we factorize the decay amplitude into soft($), hard(i/), and harder (C) dynamics 
characterized by different scales, [HI [10] 

A~ Jdx 1 dx 2 dx 3 b 1 dbib 2 db 2 b 3 db 3 Tr C{t)^ B (x 1 ,bi)^ D (x 2 ,b 2 )^ D (x 3 ,b 3 )H(x i ,bi,t)S t (x i ) e~ s( - t] . 

(4) 

In above equation, bi is the conjugate space coordinate of the transverse momentum k^r, and t is 
the largest energy scale. C is Wilson coefficient, and $ is the wave function. The last term, e~ s ^\ 
contains two kinds of contributions. One is due to the resummation of the large double logarithms 
from renormalization of ultra-violet divergence lni&, the other is from resummation of double 
logarithm In 2 b from the overlap of collinear and soft gluon corrections, which is called Sudakov 
form factor. The hard part H can be calculated perturbatively, and it is channel dependent. More 
explanation of above formula and review about PQCD can be found in many reference, such as 

[ainniin]. 

As a heavy meson, the B meson wave function is not well defined, neither is D meson. In heavy 
quark limit, we take them as: 

<S> B (x,b) = ^=[P + M B ] 75 <{ )B (x,b), (5) 

$ D (x,b) = -j=: l5 [P + M D }<l) D (x,b). (6) 

The Lorentz structure of two mesons are different because the B meson is initials state and D 
meson is final state. 

The effective Hamiltonian b — > q{q = d, s) is given by |12j : 

n cS = ^{^[d^o^/x) + c 2 Gu)o 2 c ( M )] +v ug v: b [c 1 ( f i)ouri + cmozQj) 

10 

-VZVttYrCiMOifr)}, (7) 

i=3 
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where Ci(fj,)(i = 1, • • • , 10) are Wilson coefficients at the renormalization scale fi and the four quark 
operators Oi{i = 1, • • • , 10) are 



0\ = 


(biCj)v-A(cjqi)v-A, 


oi = 


(biCi)v-A(cjqj)v-A, 


Of = 


(biUj) V -A(ujqi)v-A, 


oi = 


(biUi)v-A{ujqj)v-A, 


o 3 = 


(kqi)v-A T,qi9jQj)v-A, 




(b~iqj)v-AY,q{<lj<li)v-A, 


o 5 = 


(biq^V-A T, q i9jQj)v+A, 


o 6 = 


(biqj) V -A Y.qiQjQih+A, 


o 7 = 


§(Mi)v-A Eg e q(<ij<lj)v+A, 


8 = 


U^Qj)v-A Eg e q (qjqi)v+A, 


o 9 = 


UkQi)V-A Eg e q(Qj1j)v-A, 


Oio = 


-- Ubiqj)v-AY,q e q{Wi)V-A 



(8) 



Here i and j are SU(3) color indices; in O3 io the sum over q runs over the quark fields that are 
active at the scale \i = 0(to{,), i.e., q G {u, d, s, c, b}. For Wilson coefficients, we will also use the 
leading logarithm summation for QCD corrections, although the next-to -leading order calculation 
already exists [12]. This is the consistent way to cancel the explicit fx dependence in the theoretical 
formulae. 




Figure 2: The leading order Feynman diagrams for Bd — * D D process in PQCD approach 

According to the effective Hamiltonian in eq.(0[8|), the lowest order diagrams of B — ► D°D° are 
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drawn in Fig. [2J We first calculate the usual factorizable diagrams (a), (b), (c) and (d). For the 
(V — A)(V — A) operators, their contributions of (a) and (c) are always canceled by diagrams (b) and 
(d) respectively because of current conservation. For the {V — A)(V + A) operators, these diagrams 
can not give contribution, either. That's to say, factorizable diagrams have no contribution. For 
non- factorizable diagrams (e), (f), (g) and (h), we find the hard part for (V — A)(V — A) operators 
are same to (V — A)(V + A) operators. We group the contribution of diagrams (e) and (f), denoted 
by M a , as follows: 

UirC F Ml f 1 



M a [Ci] = —j==— J dx 1 dx 2 dx 3 I b\db\b 2 db 2 <f>B(xi,bi)4>D(x2,h)4'D{x3 1 bz) 







x j Xl + x 2 + (2x 3 -x 2 )r 2 C i {t 1 a )E(t 1 a )hi 1 \x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

+ ^-x 3 + (2x 1 -2x 2 + x 3 )r 2 ^C t (t 2 a )E(t 2 a )h^{x 1 ,x 2 ,x 3 ,b 1 ,b 2 )^, (9) 

where Cf = 4/3 is the group factor of SU(3) C gauge group, and Cj is Wilson coefficient. The 
function E m is defined as 

E(t) = a s (t) e - s B(t)SD(t)-S D (t) ; ( 1Q ) 

and Sb, Sd result from Sudakov factor and single logarithms due to the renormalization of ultra- 
violet divergence. The functions h a is the Fourier transformation of virtual quark and gluon prop- 
agators. It is defined by 

h^(xi,x 2 ,x 3 ,b 1 ,b 2 ) = 



^B^(M BV /x 2 x 3 (l-2r^b 1 )UM BV / x 2 x 3 (l-2r^ b 2 )0(h - b 2 ) 

' K (MBF aU) h), for Fl >0 

+ (61 <-» 62) \ x j _._ r —— _ _y> ) , (11) 



%T%>(M B J\Ffa\ h), for Ffa<Q 



with: 



F a 2 (1) = -x 1 -x 2 -x 3 + x 1 x 3 + x 2 x 3 + {x 2 + x 3 -xix 3 -2x 2 x 3 )r 2 ; (12) 
F a(2) = X 2^3 - x\x 3 + (xix 3 - 2x 2 x 3 )r 2 . (13) 

In above equation, Hq 1 ^^) = Jo (z) + iYq(z). In order to reduce the large logarithmic radiative 
corrections, the hard scale t in the amplitudes is selected as the largest energy scale in the hard 
part: 



t{ = max(Mfl \/\ F a(j) I ' ^3 \J (1 - 2r 2 )x 2 x 3 , 1/h, l/b 2 ). (14) 



Analogically, we can get the which comes from the contribution of diagrams (g) and (h): 



M b [Ci] = — -j==A J dxidx 2 dx 3 J bxdbxbzdbz (t> B {xiM)4>D{x2M)<t>D{xzM) 
1 - x 3 + (2 + 2xi - 2x 2 + x 3 )r 2 ] C^ 1 )^ 1 )/^ {x u x 2 , x 3 , b x ,b 2 ) 



+ 



Xl + x 2 - 1 + (-2 - x 2 + 2x 3 )r 2 ] Ci{tl)E{tl)hf \xi, x 2 , x 3 , h, 6 2 )}, (15) 
and the functions are denned as: 
h < b j \xi,x 2 ,x 3 ,b 1 ,b 2 ) = 



yH^MB^l - x 2 - x 3 + x 2 x 3 + (x 2 + x 3 - 2x 2 x 3 )r 2 b%) 

x Jo(Ms V 1 ~ x 2 - x 3 + x 2 x 3 + (x 2 + x 3 - 2x 2 x 3 )r 2 b 2 )9(h - b 2 ) 

K (M B F bU) h), 
JH«(M B .A^J6 1 ), for F 2 < 0, 



for if ( . } >0 



F, 



6(1) 
6(2) 



-1 - + x 2 x 3 + (xix 3 - 2x 2 x 3 )r 2 , 



F 2 ^ = 1 — xi — x 2 — x 3 + xix 3 + x 2 x 3 + (x 2 + x 3 - xix 3 - 2x 2 x 3 )r 2 , 



(16) 

(17) 
(18) 



t{ = max(M B ^/|F 2 (i) |, M B y/l - x 2 - x 3 + 2:3X5 + (x 2 + x 3 - 2x 2 x 3 )r 2 , 1/61, l/6 2 ). (19) 
So, the decay amplitude of decay -B^ — > D°D° can be read as: 



= VSK d M a [C 2 ] - V t * b V td M a [C B + C 7 ] + V* b V ud M b [C 2 ] 
= V^V^-VttVuPi 
= V* h V td P l (l + z 1 e i ^ +5 ^\ 



V t * b V td M b [C 5 + C 7 ] 



(20) 



where (3 is weak phase angle defined in Eq.([T]), and 5\ is the strong phase, which plays an important 
role in studying CP-violation. In above calculation, we denote that 



Ti = M a [C 2 ]-M b [C 2 ], 

Pi = M a [C 5 + C 7 ] + M b [C 5 + C 7 ] + M b [C 2 ], 



(21) 



and 



v: b Vcd 

v t * b v td 



Tj. 
Pi 



(22) 



which describes the ratio between tree diagram and penguin diagram. The corresponding charge 
conjugate decay is 



Ai = Vt b V^P L (l + z 1 e 



)• 



Therefore, the averaged decay width T for B° — ► D°D decay is then given by 



Y(B° -► D°D°) 



C 2 M 3 

-^—^ (1 - 2r 2 )| TO d Pi| 2 | 1 + ^ + 2zi cos /3 cos 5 1 \ . 

1Z07T 



(23) 



(24) 



^From this equation, we know that the averaged branching ratio is a function of CKM angle f3, if 
z\ ^ 0. Derived from Eq. ([20|) and Eq. (f23j) . the direct CP-violation can be formulated as: 



|2_U _n I 2 



-2^i sin (3 sin 5i 



(25) 



\A B ^ DQ jp\ 2 + |^4_ ^^o^J 2 1 + z\ + 2zi cos/3cosi5i 
For — > D [) D" and its conjugate decay, we write the decay amplitudes and rearrange them 



as: 



A 2 



A 2 



V; b V cs M a [C 2 ] - V t lV ts M a [C 5 + C 7 ] + V: b V us M b [C 2 ] - V t * b V ts M b [C 5 + C 7 ] 
V: b V us M b [C 2 ] - V t tV ts \M a [C 5 + C 7 ] + M b [C 5 + C 7 ] - ^^M a [C 2 }\ 

L V tb V ts > 

V: b V us T 2 - V t * b V ts P 2 

v: b v us r 2 [i + Z2 e<-^ 

V ub V*T 2 \l + z 2 e^ +5 ^ 



where T 2 , P 2 and z 2 are defined as: 



(26) 
(27) 



T 2 = M b [C 2 ], 

P2 = M a [C 5 + Cr} + M b [C 5 + C 7 ] 
T 2 



v t tv ts 



M a [C 2 



z 2 



V ub V us 



P-2 



So, the averaged decay width and direct CP violation can be formulated as: 



T(B S -> D°D°) 



C 2 M 3 o 

-^(1 - 2r 2 )\V ub V: s T 2 \ 2 (l + z 2 + 2z 2 cos5 2 cos 7 ), 



A d S r P {B s -> £> U D' 



07=^ 



B S ^D°D 



2z 2 sin 7 sin #2 



2 1 + Zo + 2z 2 cos 7 cos (52 



(28) 

(29) 
(30) 



QCD 



In our calculation, we set m c m£>, just because mo — m c Aqcd and 
quark limit. 



in the heavy 



Table 1: Amplitudes (10" 3 GeV) of B d -» D°D and B s -» D°D . 





Bd - Z^ 


Bs - L»0S° 


T(e)+T(f) 


68 + 17 i 


66 + 27i 


P(e) + P(f) 


0.80 + 0.23 t 


0.77 + 3.68 1 


T(g)+T(h) 


9.81 - 2.99 i 


14.0 - 0.6? 


P(g)+P(h) 


0.08 - 0.02 i 


-0.01 + 0.01 i 



3 Numerical Results 



For B meson, the distribution amplitude is well determined by charmless B decays [HI [TO], which 
is chosen as 

M 2 x 2 



1 



H 2 



(31) 



4>b(x, b) = Nbx 2 (1 — x) 2 exp 

with parameters oj b = 0.4 GeV, and Nb = 91.745 GeV which is the normalization constant using 
fs = 190 MeV. For B s meson, we use the same wave function according to SU(3) symmetry, where 
io b = 0.4 GeV, N Ba = 119.4 GeV and f Bs = 230 MeV. 

Since the c quark is much heavier than the u quark, the c quark shares more momentum, and 
this function should be asymmetric with respect to x = 1/2. The asymmetry is parameterized 
hy dp. Similar to the 6-dependence on the wave function of B meson, for controlling the size of 
charmed mesons, we also introduce the intrinsic ^-dependence on those of charmed mesons. Hence, 
we use the wave function of D meson as 1 13 1 



<j>D(x,b) 



f D x(l-x) 1 + ojj(1-2x) 



exp 



1 



(co D bf 



(32) 



We use an = 0.7 and ujd = 0.4 in above function. Other parameters, such as meson mass, decay 
constants, the CKM matrix elements and the lifetime of B meson are list [T| 114] : 

M B = 5.28 GeV, M Bs = 5.36 GeV, M D = 1.87 GeV, f D = 210 MeV, 
\V ud \ = 0.974, \V ub \ = 4.3 x 10~ 3 , \V cd \ = 0.23, \V cb \ = 41.6 x 10" 3 



\Vtd\ = 7.4 x HT 3 , \V tb \ = 1.0, |K 
WtJ = 41.6 x 10~ 3 



0.226, \Vcs\ = 0.957, 



t b o = 1.54 x 10" 12 s, t b o = 1.46 x 10" 12 s. 



(33) 



With these parameters fixed, we calculate the decay amplitudes of the B u — > D°D and B 



D°D° decays in Table [TJ From the table, we notice that the main contribution comes from the tree 




Figure 3: The branching ratio of B s — > D D changes with CKM angle 7. 



diagram (e) and (f). And our predictions for the branching ratio of each mode corresponding to 
(3 = 23° and 7 = 63° are listed, 

BR(B d D°lf) = 2.3 x 1(T 5 ; 

BR(B S -» D°D°) = 6.8 x 10" 4 . (34) 

In Fig. El we plot the branching ratio of B s — ► D°D° with different 7. In this figure, we find the 
branching ratio is not sensitive to CKM angle 7. For the experimental side, there are only upper 
limits given at 90% confidence level for decay B d — > D^T^ , 

5i2(.B d D !} ) < 6.0 x 10 -5 ; SarSarpS] 

-» L>°D°) < 4.2 x 10~ 5 . 5eZZe[IS] (35) 

Obviously, our result is consistent with the data. For Bd — ► D°D^ decay mode, zi is about 6.5, 
and the strong phase S\ is 34°, so A^p is about —6% with the definition in Eq. (|25|) . As decay 
mode B s — > D°D° is concerned, ^2 is about 205 and ^2 = 155°, and the relation between direct 
CP violation and 7 is shown in FigUl From the figure, we read the CP asymmetry is about 0.4%, 
which is rather tiny. It is necessary to state that the z\ and Z2 are not the true ratio between tree 
contribution and penguin, because mathematical technics are used in Eq. (|20p and (|27p . 

In addition to the perturbative annihilation contributions, there is also a hadronic picture for 
the B^ — > D°T^ , named soft final states inter action [T7]. The B meson decays into D + and D~ , 

9 
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Figure 4: The direct CP-violation of B s — ► D°D changes with CKM angle 7. 

the secondary particles then exchanging a p meson, then scatter into D°D° through final state 
interaction afterwards. For B s decay, the B s meson decays into D^t and D + then scatters into 
D°D° by exchanging a Kaon. But this picture cannot be calculated accurately because of lack of 
many effective vertexes, and we will ignore this contribution here, though it may be important [17J. 

There are many uncertainties in our calculation such as higher order corrections, the parameters 
listed in Eq. (|33p and the distribution amplitudes of heavy mesons. We will not discuss uncertainty 
taken by high order correction as we only roughly estimate the branching ratios and CP asymme- 
tries, though high order corrections have been done for some special channels |18t [T9] and showed 
15 — 20% uncertainty. The parameters in Eq. (|33p . fixed by experiments, are proportional to the 
amplitudes, so we will not analyze this kind uncertainty either. In our calculation, we find that 
the results are sensitive to the distribution amplitudes, especially to that of D meson. Since the 
heavy D wave function is less constrained, we set au € (0.6 — 0.8) GeV and ud G (0.35 — 0.45) GeV 
to exploit the uncertainty. Table [2] shows the sensitivity of the branching ratios to change of u^, 
ojb and ao- It is found that uncertainty of the predictions on PQCD is mainly due to w^i, which 
describes the behavior in end-point region of D meson, however it is very hard to be determined. 
Considering the experimental upper limit, our results favor large (Jb, large ujd and small ao- 

At last, we give the prediction of branching ratios with err bar as follows: 

B^anf) = (3.8±°;ii;;ltg:l) * ( ^v^v^ooev ) 2 ; 
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Table 2: The sensitivity of the decay branching ratios and CP asymmetries to change of Ub, lod 
and an 





BR(B d -> D°D°) 
(xl0~ 5 ) 


BR(B S -> D°D°) 
(xlO- 4 ) 


Ac r P {B d -> D°D°) 
(%) 


A*%(B S -> 
(%) 




uj h {B \ B s ) 
0.35 \ 0.45 
0.40 \ 0.50 
0.45 \ 0.55 


4.3 
3.8 
3.2 


7.8 
6.8 
5.9 


-7.2 
-5.3 
-5.8 


0.4 
0.4 
0.4 




0.35 
0.40 
0.45 


5.0 
3.8 
2.2 


9.7 
6.8 
4.2 


-4.2 
-5.3 
-7.8 


0.3 
0.4 
0.5 




a D 
0.6 
0.7 
0.8 


3.2 
3.8 
4.3 


5.9 
6.8 
7.8 


-6.9 
-5.3 
-6.1 


0.4 

0.4 
0.4 





We believe that the Bj — > Z) Z) will be measured soon because this ratio is just below the upper 
limit, and Bj — ► D°D° will be measured in LHC-b in next year as a channel to cross check the 7 
measurements. 

4 Summary 

With heavy quark limit and hierarchy approximation Xqcd *C mp <C ms, we analyze the 
B -> Z)°Z) and B s -> Z)°S decays, which occur purely via annihilation type diagrams. As a 
roughly estimation, we calculate the branching ratios and CP asymmetries in PQCD approach. 
The branching ratios are still sizable. The branching ratio of B —> D°D° is about 3.8 x 10~ 5 , which 
is just below the experimental upper limited result [T5J [16], and we think that it will be measured 
in near future. For B s — * D°D° , the branching ratio is about 6.8 x 10 -4 , which could be measured 
in LHC-b. From the calculation, it is found that this branching ratio is not sensitive to angle 7. 
In these two decays, there exist CP asymmetries because of interference between weak and strong 
interaction, though they are very small. 
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